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We have studied magnetic Feshbach resonances in an ultracold sample of Na prepared in the 
absolute hyperfine ground state. We report on the observation of three s-, eight d-, and three g- 
wave Feshbach resonances, including a more precise determination of two known s-wave resonances, 
and one s-wave resonance at a magnetic field exceeding 200 mT. Using a coupled-channels calculation 
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leading to the Feshbach resonances in terms of the asymptotic-bound-state model is presented. 
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I. INTRODUCTION 

The observation of Feshbach resonances in ultracold 
atomic gases has opened a wealth of exciting ex- 

periments, as their presence allows tunability of the two- 
body interaction strength and coherent association of ul- 
tracold molecules. Feshbach resonances are caused by 
the coupling between an atom pair and a molecular state 
in a closed channel potential [5|. The most commonly 
used ones are magnetically induced resonances that orig- 
inate from a difference in the magnetic moment of the 
molecular state and the atomic threshold Q. The loca- 
tion of the resonances is very sensitive to the interaction 
potentials. Much effort has gone into constructing ac- 
curate two-body potentials in order to predict Feshbach 
resonances and scattering properties. Conversely, precise 
determination of Feshbach resonances can be regarded as 
a sensitive spectroscopic tool to map out the molecular 
spectrum of the least bound states. 

For Na only three Feshbach resonances have been re- 
ported so far, namely two for an ensemble prepared in the 
absolute hyperfine ground state \F,mp) = |1,1) [l[ and 
one in the hyperfine substate |1,— 1) H- AH three are 
s-wave resonances, i. e. caused by molecular states with 
rotational quantum number 1 = 0. In parallel, Raman 
spectroscopy on a molecular beam of sodium dimers to 
locate the least-bound rovibrational states of the ground 
state potentials X^T,'^ and a'^E+ has been performed 0- 
[Toj . Scattering lengths @ as well as Feshbach resonance 
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positions [T3l were extracted. In addition, Refs. [TTl, 113 
presented an extensive study of the lowest triplet po- 
tential a"^S^ by two-color photoassociation spectroscopy 
on magneto-optically trapped Na atoms. However, more 
precise Feshbach spectroscopy of Na is needed to improve 
the knowledge of the interaction properties. 

Here, we report on the observation of seven d-wave and 
three g-wave Feshbach resonances (i. e. caused by molec- 
ular states with I = 2 and I = 4, respectively) in the 
magnetic field range between 49 mT and 91 mT, and a 
more precise determination of two s-wave Feshbach res- 
onance positions first observed in Ref. [l[. In addition, 
we have observed an s- and d-wave Feshbach resonance 
at 205.4 mT and 159.0 mT, respectively. Using coupled- 
channels calculations we have improved the Na2 ground- 
states potentials by combining our new Feshbach reso- 
nance data with the known spectroscopy data. We have 
also applied the recently developed asymptotic-bound- 
state model (ABM) [Ullll for the relevant / = 0, 2 and 
4 molecular states of Na2, providing a simple description 
of the data. Note that g-wave Feshbach resonances have 
so far only been observed for Cs [l^ [IB] , for which even 
molecular states with I = 8 have been populated [13, [13 ■ 

This paper is organized as follows. In Sec.|ll]we out- 
line our experimental procedure for the preparation of 
a Na Bose-Einstein condensate (Sec. Ill A|) and the Fesh- 
bach spectroscopy measurements (Sec. lIIB) ). In Sec.lIIII 
the theoretical model for the near threshold molecular 
levels of Na2 is presented, introducing the interaction 
Hamiltonian (Scc. lIII A|) . followed by a short descrip- 
tion of the ABM and its application to the Feshbach 
resonance data (Sec lHIBj) . We present results of the 
coupled-channel calculation deriving improved ground- 
state potentials of Na2 (Sec. lIII C|) . We summarize, give 
the scattering lengths obtained from the new potentials. 
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FIG. 1: (color online) Feshbach spectroscopy measurements for Na prepared in the lowest hyperfine substate |_F, jtif) = |1, 1), 
showing the number of atoms observed after a certain hold time as function of the magnetic field B. The field range for each 
panel differs and longer hold times indicate weaker and narrower resonances. Blue, red and green colors corresponds to s-, d-, 
and g-wave resonances, respectively. Each data point is an average over 3 to 4 experimental runs; the error bars represent the 
one standard deviation statistical uncertainty. The solid lines are Gaussian fits, used to determine the loss resonance position. 
Note that for panel (/) the right side of the resonance is not taken into account in the fit. 



and conclude in ScclTV 



II. EXPERIMENT 
A. Preparation of a Na Bose-Einstein condensate 

The observation of magnetically induced Feshbach res- 
onances requires the preparation of a sufficiently high 
density and sufficiently cold atomic sample, on virhich an 
external homogeneous magnetic field can be applied. We 
use a Na Bose-Einstein condensate (BEC) in an optical 
dipole trap (ODT). Our experimental procedure is based 
on evaporative cooling of Na in the |2, 2) state in a cylin- 
drical symmetric cloverleaf lofFe-Pritchard magnetic trap 

We load about 10'^ atoms from a dark-spot magneto- 
optical trap into the magnetic trap, and then spin- 
polarize at a bias offset field of 2.0 mT into the |2,2) 
state, in analogy to the strategy of Ref . [20| . Resid- 
ual atoms in |2, 1) are removed state-selectively by mi- 
crowave radiation, using an offset field of 5.3 mT. Then, 
we apply forced evaporative cooling driving the |2,2) 
to |1,1) transition for about 30 s. We first use a com- 
pressed magnetic trap with axial and radial trap frequen- 
cies (faxji'rad) = (33, 400) Hz, while during the last few 
seconds of the evaporation stage we expand the mag- 
netic trap to trap frequencies (I'ax, t'rad) = (33, 130) Hz 
in order to reduce losses from three-body recombination. 
After evaporation, we are left with 1 x 10'' atoms at a 
temperature of 2 /xK. 

At this point we load the atoms into an ODT, formed 
by two crossed laser beams with waists of 60 fim each and 



powers of 1 W and 0.5 W, respectively, and derived from 
a 1064 nm Nd:YAG laser. The resulting trap has a depth 
of about 8 /iK and frequencies (j^ax, i^rad) = (130, 310) Hz. 
By our loading procedure and subsequent rethermaliza- 
tion by elastic collisions, we routinely obtain a BEC of 
5 X 10^ atoms in |2,2). Subsequently, we transfer all 
atoms to the state |1, 1) by rapid adiabatic passage, ap- 
plying a 50 ms magnetic field sweep of 24 fiT around an 
offset field of 0.1 mT at a fixed microwave frequency. The 
typical lifetime of the |1,1) BEC with peak density of 
3 X 10^"* cm~^ in the ODT is about 5 s, limited by three- 
body recombination. 



B. Feshbach spectroscopy 

A common way to probe Feshbach resonances is to 
measure the rate of atom loss as a function of magnetic 
field B. The dynamics in a Bose gas is controlled by 
collisions that are characterized by the s-wave scattering 
length a. Near a Feshbach resonance a can be described 
as 
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aiB) = flbg 1 



B-Br 



(1) 



where Obg is the background scattering length, Bq the res- 
onance field and A the magnetic width of the resonance. 
The presence of a Feshbach resonance leads to a strong 
enhancement of loss because the dominant loss process, 
namely three-body recombination, scales strongly with a 

To locate a resonance, we ramp to a magnetic field B, 
wait for a certain hold time, after which we ramp the 
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field back to zero and use absorption imaging to mea- 
sure the remaining number of atoms. The ramp times 
are always short compared to the hold times. A com- 
pilation of our measurements is shown in Fig.[T] and the 
experimental results are summarized in TableU] We have 
observed 12 loss resonances between 49 mT and 91 mT, 
which we assign as seven d-wave and three y-wavc reso- 
nances. The s-wave resonances at 85.1 mT and 90.5 mT 
have been measured with higher accuracy than reported 
previously [l|. Guided by initial theoretical models we 
performed a restricted search at higher magnetic fields 
and located two additional Feshbach resonances, namely 
one c?-wave resonance at 159.0 mT and one s-wave reso- 
nance at 205.4 mT. 

The required hold times range from 25 ms for the 
broadest s-wave resonances up to 4 s for the very narrow 
d- and 5- wave resonances. Except for the s-wave res- 
onances, the typical width of the loss features is about 
10 /iT to 20 /iT, mostly determined by our magnetic field 
stability. We take the magnetic field value of maximal 
loss as the resonance position B'^^ . The magnetic field 
is calibrated close to each Feshbach resonance by means 
of RF spectroscopy on the 1 1 , 1 ) 1 1 , 0) transition of Na. 
The asymmetry observed in the s-wave loss features in 
Fig.[TJk) and[TJ/) is caused by molecule association when 
ramping down the magnetic field for imaging. 

The measured position of the broadest s-wave reso- 
nance of 90.51(4) mT is in agreement with the previ- 
ously reported value of 90.7(2.0) mT [H. We find the 
spacing between the lowest two s-wave resonances to be 
5.41(5) mT, also in agreement with the earlier value of 
5.4(1) mT [l|. We note that preliminary data on part of 
the Na d-wave resonances have been reported in Ref. [131 . 

III. THEORY 
A. Hamiltonian 

Feshbach resonances occur when molecular states are 
resonant with two scattering atoms. For sodium atoms 
in the ground state the molecular states are re- 

lated to the least-bound rovibrational levels of the sin- 
glet A^E+ and triplet a'^E+ Born-Oppenheimer poten- 
tials. Their rovibrational levels are labeled by the vi- 
brational and rotational quantum numbers v and I as 
well as the space-fixed projection mi of the orbital angu- 
lar momentum along the magnetic field direction. The 
basis |(t) = \SMsIMi) describes the remainder of the 
molecular wavefunction. The spins S and / are the total 
electron and nuclear spin of the molecule, respectively. 
In fact, for a pair of atoms A and B we have S = sa + sb 
and I — I A + IB, where sa/b and ia/b are the electron 
and nuclear spin of each of the atoms. Projections are 
again defined along the magnetic field direction. 

There is significant hypcrfinc interaction in Na and it 
is also useful to define the coupled basis \{SI)f,mf), in 
which the total electronic and nuclear spins are coupled 



TABLE I: Overview of the experimentally and theoretically 
obtained results on Feshbach resonances in Na prepared in 
the lowest hyperfine substate ttif) = Experimen- 
tally the position of maximum loss Bg"'' is given, which is de- 
termined by Gaussian fits to the loss features. The errors re- 
flect the one standard deviation statistical uncertainty in the 
magnetic field calibration, determined by the FWHM of the 
RF calibration spectra and from the profile fit. The middle 
columns shows the results for the Moerdijk and ABM models 
in terms of Be (see Sec. lIII B)) . The last two columns show the 
Feshbach resonance position Bq and width A from coupled- 
channels (CC) calculation, based on refined Na2 ground-state 
potentials (see Sec. lIII C|) . The brackets in the last column 
give the exponent to the power ten. The quantum numbers 
I, f, ruf characterize the bound state and are discussed in the 
text. Note that for the ?=4 states the m/ quantum numbers 
could not be assigned (see Sec. lIII C)) 
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to / = S + I, and m/ is its projection. We will la- 
bel the molecular states with quantum numbers I and 
nif, common to both bases, and if appropriate with / 
or Ms- For homonuclear dimers symmetrization ensures 
that only states with I + S + I even exist [2^ . 

In order to accurately and quantitatively describe both 
bound states and scattering processes, which include 
Feshbach resonances, we require the Hamiltonian for two 
sodium atoms in a magnetic field B (see e.g. [13, [28j). 
Here, we use an extension that allows for a second-order 
spin-orbit interaction [29j and a hyperfine-contact inter- 
action that depends on the internuclear separation R. 
The effect of this second correction was recently studied 
for Rb2 [131 • Finally, we have 

+VssiR), (2) 

where T = — fi,^V^/(2^) is the relative kinetic energy 
operator with reduced mass ^, 

a=A,B 
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is the i?-dcpendcnt hypcrfine-contact interaction, and 

Hz = ^ [gsaS^a + 9iaiza)lJ-BB/h 

a=A,B 

is the magnetic Zeeman interaction. In the limit i? — > oo, 
the i?-dependent functions aa{R) approach the atomic 
values flhf . The constants gson gia, and fis are the sodium 
electron gyromagnetic ratio, sodium nuclear gyromag- 
netic ratio, and the Bohr magneton, respectively. All 
atomic constants are taken from Ref. [3l[ except for the 
atomic mass of sodium, which is obtained from Ref. [s^I ■ 

The last two terms on the first line of Eq. [2] describe the 
Hamilton operator for the Born-Oppenheimer potentials 
Ux{R) and Ua{R)- The operators Px and Pa project on 
spin basis functions \SMsIMi) with total electron spin 
5 = and 1, respectively. Consequently, Px + Pa = 1- 
Finally, Vss{R) is a weak but non- negligible effective 
electron-spin electron-spin interaction that includes the 
second-order spin-orbit interaction and depends on the 
orientation of the internuclear axis R as well as the inter- 
nuclear separation. Therefore, Vss{R) is the only term 
in the Hamiltonian that mixes or couples partial waves 
I. It will be defined more precisely in the Sec. lIII Cl 

The Hamiltonian conserves the total angular momen- 
tum / + / at zero magnetic field, and always conserves 
its projection M = nif + mi. Neglecting the potential 

Vss{R), the Hamiltonian conserves / for B = 0, while 
for finite field only the projection m/ remains a good 
quantum number. For large magnetic field, where Hz 
is larger than Hhf{R), the total electron spin S and its 
projection Ms are approximately conserved. For sodium 
this corresponds to fields B > 30 mT. 

We focus on ultra-cold collisions, with temperatures 
well below the p-wave centrifugal barrier. Hence, only 
s-wave collisions need to be considered. This also means 
that the d- and g-wavc resonances can only be induced 
by the weak spin-spin interaction Vss{R), which allows 
coupling between the s-wave continuum and / = 2 and 
4 molecular states, respectively. Moreover, for Na atoms 
in the |1, 1} hyperfine state, = mj^ -f =2 and, 
therefore, we can focus on basis functions with M = 2. 
The spin basis states \a) for I = and 2 are listed in 
Tableim Also shown is an equivalent list in the coupled 
\{SI)f,mf) basis. With the constraint AI = 2, there 
are five / = basis states, where only one has singlet 
character, whereas for I ~ 2 there are four singlet and 18 
triplet basis states. 



B. Asymptotic-bound-state model 

In general, it is difficult to predict the location of Fesh- 
bach resonances, because of their sensitivity to details of 
the molecular potentials. This has lead to the develop- 
ment of models that require minimal knowledge of the 
molecular potentials (see e. g. [33l - l36j ). at the expense of 



TABLE 11: List of / = and 2 states for each vibrational 
singlet {S = 0) and triplet (S = 1) state, for which A4 = 
nif + mi =2 holds. The last two columns show the quantum 
numbers of the coupled basis, which holds for low magnetic 
fields. 
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accuracy and applicability. A powerful, yet computa- 
tionally light description of the near threshold molecular 
sgectrum is the asymptotic bound-state model (ABM) 
ImM- It builds on an earlier model by Moerdijk et 
al. [21| for homonuclear systems, which was extended 
by Stan et al. [13 for heteronuclear systems. By using 
binding energies and wavefunction overlaps between sin- 
glet and triplet rovibrational levels as fit parameters it 
circumvents the need of any knowledge of the molecular 
potentials while allowing for a high degree of accuracy 
once the position of a few Feshbach resonances is known. 
This type of modeling is closely related to the well de- 
veloped deperturbation theories of molecular spectra, see 
e.g. K2 nil. 

In the ABM model we neglect the R dependence of 
the hyperfine constants and set Vss{R) to zero in the 
Hamiltonian of Eq. [5J We then define the molecular 
wavefunctions \vl,a) = |'!/'^'')|cr), where IV'^'O describe 
vibrational and rotational wavefunctions of the 5 = 
singlet or 5 = 1 triplet potential, and formally construct 
the Hamiltonian matrix % with matrix elements 

{v'l',a'\H\vl,a) 

for a limited set of near-threshold bound states of the two 
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Born-Oppenhcimcr potentials and spin states \a) consis- 
tent with total projection quantum number M . 
For our homonuclear system we find 

T~i-v' V a' ,vla ~ £ y ^v v' ^IV 5<j<j' (3) 

a 

where 6ij is the Kronecker delta, e^'' is the energy of the 
rovibrational level {v, I) of the singlet or triplet potential, 

and T]^^, (/) = ^^'^ overlap integral of two 

rovibrational wavefunctions. Clearly, for S = S' we have 
Vvv'iO — ^vv'- The operator Sq-«*q does not commute 
with the total electron spin 5*. In fact, it is the only part 
of H that couples singlet and triplet levels. 

The eigenvalues of the Hamiltonian TL arc the molec- 
ular bound states. We can label these levels by M, trif, 
and / as well as their B = quantum number /. As 
function of magnetic field some of these molecular bound 
states cross the field-dependent |1, 1) -I- |1, 1} dissociation 
limit. Such a crossing position is denoted by Be, which 
principally differs from the Feshbach resonance position 
Bq (as defined in Eq.[T]) because finite coupling between 
the molecular bound states and the scattering channel 
shifts the resonance position Q . For the goal of the ABM 
of giving a simplified picture the shifts are not important 
in the case of Na2. We adjust the energies e^'' and over- 
laps ry^f / (/) to fit the molecular spectrum to the observed 
Feshbach resonances. 



1. Moerdijk model 

We first neglect singlet-triplet coupling by applying the 
following approximation in Eq. [3j 

a 

such that the hyperfine interaction commutes with the 
total electron spin S. We refer to this approximation 
of ABM as the "Moerdijk model" [2l|. In the Moerdijk 
model the only fit parameters are the energies e^'', as the 
overlap parameters r]^^, (/) do not appear. 

We find that all the observed Feshbach resonances can 
be described by a single vibrational triplet level below 
the |1, 1) + |1, 1) dissociation threshold. In fact it is the 
Wq = 14 level with I = 0, 2 and 4 From a least- 

squares fit including all resonances we obtain for the ro- 
tational structure of this vibrational level three fitted en- 
ergies: elf/h = -4976 MHz, el'^/h = -3679 MHz and 
e\f/h = —765 MHz, where the energies are given with 
respect to the |1, 1) + |1, 1} dissociation threshold. The re- 
sulting molecular spectrum is shown in Fig. [2] and demon- 
strates that the Moerdijk model describes the spectrum 
of Feshbach resonances already quite well. The corre- 
sponding Feshbach resonance positions are given in Ta- 
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FIG. 2: (color online) Moerdijk model of the Va ~ 14 triplet 
vibrational level with M — 2, and I = (thick blue lines), 
/ — 2 (thin red lines) and I = 4 (dotted green lines) rota- 
tional states. The dashed black line is the atomic |1, 1) + 11, 1) 
threshold. The observed Feshbach resonances are indicated 
by the filled squares (s-wave resonances), circles (d-wave res- 
onances) and triangles (g-wave resonances). Note that the 
three g-wave resonances overlap with one d-wave resonance 
in this graph. On the left side of graph the 5 = quantum 
number / is indicated for the different I states, whereas inside 
the graph the approximate quantum number Ms is denoted. 

bleU under "Moerdijk" . The deviations with the experi- 
mental values are at most 0.3 mT, much smaller than the 
spacing between the resonances, and allows for an un- 
ambiguous labeling of the molecular states causing the 
resonances. The observation of three g-wave Feshbach 
resonances near 51 mT close to the d-resonances can only 
be explained by the presence of / = 3 {Ms = 0) states. 
However, an TO/-assignment is ambiguous as singlet- 
triplet coupling and spin-spin interaction plays a signifi- 
cant role for these nearly degenerate states. Thus no to/ 
label is given in Table|T]for these cases (see Sec. lIII Cl for 
further discussion). 

In Fig. [5] the level structure is continued beyond the 
|1, 1) -I- |1, 1) dissociation threshold into the continuum. 
All these levels have finite life times. The group of Z = 4 
levels, which for all B have an energy above the dissoci- 
ation limit, could give rise to shape resonances. 



2. ABM 

To improve the description of the Feshbach data one 
has to include the singlet-triplet coupling and consider all 
rovibrational levels of the singlet and triplet potentials 
with binding energies less than a few times the Zeeman 
and hyperfine splitting. For Na2 these are the Va = 14, 
15 and vx = 64, 65 rovibrational levels [9l. [lOl. [l^ . 

There are eight fit parameters (four energies and four 
overlap integrals) for / = 0, for which there are only three 
resonances to constrain their values. If we assume that 
the energies and overlaps of different I states arc indepen- 
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dent, we have in total 32 fit parameters to fit only 14 reso- 
nances. Consequently these fits are underconstrained. In 
principle, however the number of fit parameters can be 
reduced by noting that binding energies and overlaps for 
I > follow from that for ^ = by u sing model potentials 
or the accumulated phase method However, as 

we want to use the ABM model in its simplest form we fit 
the £^'' as well as the rj^^, (l) parameters independently. 

We can improve upon the fit based on the Moerdijk 
model. We have identified a class of states that arc 
mostly insensitive to singlet-triplet mixing. The loca- 
tions of the resonances caused by / = 4 states only 
depend to well within the experimental uncertainty on 
the parameters. Thus, we can improve the level en- 
ergies compared to those obtained within the Moerdijk 
model, as we can exclude resonances that are sensitive to 
singlet-triplet mixing. We find slf/h = -4991(1) MHz 
and e'{^/h = -3682(2) MHz. The Feshbach resonance 
positions caused by these f = 4: states are given in Ta- 
bleUunder "ABM" . The calculated positions of the / = 4 
d-wave resonances deviate only by 0.05 mT or less from 
the experimental values. 

To apply the ABM model to Feshbach resonances that 
are sensitive to the vx = 64 and 65 singlet states, and in- 
directly to the Va = 15 triplet state, we include previous 
experimental data to reduce the number of fit param- 
eters. First, the e\'^ parameter can be extracted from 
photoassociation spectroscopy measurements [12| using a 
state that is insensitive to singlet-triplet mixing. One ob- 
tains el'°/h = -^2014(10) MHz. From magnetic field de- 
pendent molecular beam spectroscopy experiments on an 
avoided crossing between the Va = 14 (/ = 2, Ms = —1) 
and Vx = 64 (/ = 2, Ms = 0) states around B = 170 mT 
we extract el'"/h « -11 GHz and 77^4,14(0) ?» 0.85. 

We then find a satisfactory fit for all s-wave resonances 
with the values 0.47, 0.79 and 1.0 for the remaining three 
overlap parameters '764,15(0), '765,14(0) and ??65,i5(0). The 
fitted s-wave resonance positions are given in TablelH 
The fit is not unique, highlighting the difficulty of us- 
ing a deperturbative ABM approach for this system. 

For the d-wave resonances we follow a similar ap- 
proach, extracting 7?64,i4 (2) « 0.85 from Ref. [13], and 



using the rotational splitting from Ref. 



h X 1556(38) MHz, 



1,2 

-15 



-1,0 
-15 



.0,2 
-65 



.0,0 
-65 



-0,2 
-64 



-0,0 
-64 



h X 504(38) MHz and 



h X 465(34) MHz. With the three remain- 
ing d-wave resonances we obtain rj^l 15(2) w 0.44 and 
'765 14(2) ~ 0.47 and calculate the c?-wave resonances po- 
sition (see Table[I]). 

Finally, for the three g- waves resonances, which are 
caused by a group of Z = 4, / = 3 {Ms = 0) states, we 
only find a weak dependence on the vx = 64, 65 and 



.0,0 

-64 



= 15 levels. We use the rotational splitting 



0,4 
64 



5064(45) MHz from Ref. [9|, and find the I = A 
level of the vx = 65 and Va = 15 levels assuming the 
ratio between the / = 2 — / = and 1 = 4 — 1 = 
rotational splitting of the vx = 64 level. From a least- 
square fit we obtain the energy parameter of the Va = 14 



triplet state, e^l^ /h = —777(2) MHz, and calculated the 
g-wave resonance positions (see Table|T]). 

In summary, we can extract within the ABM the bind- 
ing energies of the Va = 14 level by selecting those reso- 
nances that are not sensitive to singlet-triplet coupling. 
To obtain the other resonances we need to include spec- 
troscopy data to restrict the number of fit parameters. 
However, the obtained set of fit parameters is not unique. 
For instance, wc have found that rj^^ and Egg are cor- 
related, thus choosing a different 7755 ic,(l) would imme- 
diately result in a different ^65. In fact, from the I = 
fa = 15 and vx = 65 wavefunctions obtained by the 
coupled-channel calculation (sec Sec lIII C]) we find that 
77g5 ]^5(0)=0.60, which significantly deviates from our fit. 

C. Coupled-channels calculation 

The description based on the Moerdijk and ABM mod- 
els agrees with the experimental resonance locations to 
well within one percent. However, the ABM is only valid 
over the limited range of energies spanned by the included 
vibrational levels and can not be used directly without 
further assumptions to find scattering lengths. Any ex- 
trapolation out of this limited energy region is not reli- 
able. Additionally, the approach needs assumptions on 
other parameters like overlap integrals to obtain a num- 
ber of free parameters that is at least equal to the number 
of observations. Thus the overall validity of this approach 
should be checked by a more general theoretical model. 

Wc use a coupled-channels approach based on the 
Hamiltonian in Eq.[2] with i?-dcpendent interaction po- 
tentials and coupling terms. The i?-dependent functions 
are found in a fit to all known Feshbach resonances and 
existing spectroscopic data of more deeply bound rovi- 
brational levels of the AT^S^ and a'^Sj potentials 
[1^ . Data from conventional spectroscopy by Kush and 
Hessel [s^ and by Barrow et al. [13] for the ground state 
A^E+ potential, and by Li Li et al. [4l| for the lowest 
triplet state a'^S+ potential are also included. The fits 
based on the coupled-channels solutions will then lead to 
a consistent description of the sodium dimcr that is valid 
over the full depth of the A^I]+ and a'^E+ potentials 
and will be a well justified system to predict scattering 
properties or dynamics like spin exchange. 

The hyperfine-contact interaction and Vss{R) in Eq.[2] 
require further discussion. The hyperfine interaction 
aa{R) is R dependent and accounts for the electronic 
distortions of one atom by the other. Wc use the simple 
Ansatz 



aa{R) 



1 



g(fl-flo)/Ai? _^ I 



(4) 



which represents a switching function around Rq with 
width AR to go from the atomic value Ua^hf at large 
R to the molecular value aQ,hf (1 + c/) at small R. For 
our homonuclear diatom we have aA{R) = clb{R)- The 
parameter values are determined by the fitting process. 
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The effective spin-spin interaction, which couples the 
partial waves I, is given by 



(5) 



where the operator Sz is the total electron spin projec- 
tion operator along the internuclear axis Z and 



A(i?) 



1 

i?3 



asoe 



-bso-R 



(6) 



The function A(i?) consists of two terms. The first rep- 
resents the magnetic dipole-dipole interaction, the sec- 
ond is modeling the second order spin-orbit contribution. 
When A(i?) and R are given in atomic units, a is the fine 
structure constant. For sodium this Ansatz was first used 
in Ref. [l^l, where they found that a non- negligible aso 
is needed to describe the photoassociation data. 

We represent the X^'E'^ and a'^Sj Born-Oppenheimer 
potentials using a separate functional form for three re- 
gions of the internuclear separation: a short range repul- 
sive wall for R < R-mn, an asymptotic long range region 
for R > i?out, and an intermediate deeply bound region 
in between. The analytic form of the potentials in the in- 
termediate region, Uib.{R), is the finite power expansion 



(7) 



in the nonlinear i?-dcpendent function ^(i?) given by 

R ^ Rm 



R + bR„ 



(8) 



Here, the upper limit n and the arc fit parameters and 
we choose b and i?™ such that only a small number of 
are needed. The separation R^ is chosen near the value 
of the equilibrium separation. For R < Rmn the potential 
is extrapolated with 



Usr{R) =ui+ U2/R' 



(9) 



where ui and U2 are adjusted to create a continuous and 
differentiable transition at i?inn. We use s = 6 for both 
X^S]+ and a'^S+ potentials. 

For R > i?out we adopt the long range form 



Ulr{R) = -Ce/R''-Cs/R''-Cw/R'''TE,^ci.{R) , (10) 

where the exchange contribution changes the sign be- 
tween singlet and triplet state and the form is given by 

m 



£^oxch(i?) = Aexi?^exp(-/3i?) 



(11) 



We start the optimization process of the potentials 
with functions known from earlier work [l0| in iteration 
loops containing two steps. First, the location of the 
Feshbach resonances and the photoassociation data from 
Ref. [ll| arc used to fit the long range behavior of Eq.lTUl 



TABLE III: Parameters of the analytic representation of the 
X^S^ state potential. The energy reference is the dissocia- 
tion asymptote. Parameters with an asterisk * ensure smooth 
continuous extrapolation of the potential at iiinn. 



R < 



: 2.181 A 



«1* 




-0.785318644 x 10"' cm" 


-1 


U2* 




0.842586535 x 10*^ cm"!^ 


16 




^inn 


<R< Rout = 11.00 A 






b 






-0.140 


Rm 






3.0788576 


A 


ao 




-6022.04193 cm" 


-1 


a\ 




-0.200727603516760356 x 10^ cm" 


-1 


0,2 




0.302710123527149044 x 10^ cm" 


-1 


13 




0.952678499004718833 x 10"' cm 


-1 


a4 




-0.263132712461278206 x 10^ cm 


-1 


fis 




-0.414199125447689439 x 10^ cm 


-1 


aQ 




0.100454724828577862 x lO'' cm^ 


-1 


ai 




0.950433282843468915 x 10^ cm" 


-1 


ag 




-0.502202855817934591 x lO'^ cm" 


-1 


ag 




-0.112315449566019326 x lO'^ cm" 


-1 






0.105565865633448541 x 10^ cm" 


-1 


ail 




-0.626929930064849034 x 10* cm^ 


-1 






-0.134149332172454119 


x 10^° cm 


-1 


ai3 




0.182316049840707183 


x 10^° cm 


-1 


0-14 




0.101425117010240822 


X lO-'-' cm" 


-1 


0-15 




-0.220493424364290123 


x 10^^ cm 


-1 


016 




-0.406817871927934494 


x 10^^ cm" 


-1 


017 




0.144231985783280396 


x 10^2 cm" 


-1 


018 




0.379491474653734665 


x 10" cm" 


-1 


019 




-0.514523137448139771 


X 10^2 cm" 


-1 


020 




0.342211848747264038 


X 10^2 cm 


-1 


021 




0.839583017514805054 


X 10^2 cm 


-1 


022 




-0.131052566070353687 


X lO^'^ cm" 


-1 


023 




-0.385189954553600769 


X 10^^ cm" 


-1 


024 




0.135760501276292969 


X lO^'^ cm" 


-1 


025 




-0. 108790546442390417 


X 10^'^ cm" 


-1 


026 




0.282033835345282288 


X 10^2 cm" 


-1 


RoMt < R 








0.0 cm" 


-i 


c% 




0.75186131 xlO'' cm-i^ 




Cg 




0.1686430x10'' cra'^J 


18 


Clo 




0.3081961 xlO^o cm-iA 


10 


Acx 




0.40485835x10^ cm-^A" 


-7 


7 






4.59105 


fi 






2.36594 A" 


-1 


Derived constants: 


equilibrium distance: 


Rf= 3.07895(5) A 




electronic term 


energy: 


T^=-D^= -6022.0420(40) cm-^ 





With the new potential functions from this first step all 
pure singlet and triplet rovibrational levels related to 
the Feshbach and photoassociation data are calculated 
and the results are added to the spectroscopic data of 
Refs. [Tol.l39l-l4ll|. Then the second fitting step yields the 
improved potentials including all regions in R. With this 
result the first fit was repeated from which again better 
rovibrational levels will be derived for the overall poten- 
tial fit. This iteration was continued until a self consis- 
tent status was reached. Because the g-wave resonances 
require long computation times, they are not included in 
the fit, but in each iteration loop the consistency of these 
data was checked. 

In the fit the position of the Feshbach resonances are 
determined by the maximum of the clastic collision rate 
at a collision energy between 1 and 2 /iK, correspond- 
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TABLE IV: Parameters of the analytic representation of the 
potential curve of the a^Ej state. The energy reference is the 
dissociation asymptote. Parameters with an asterisk * ensure 
smooth continuous extrapolation of the potential at i?inn- 



: 4.2780 A 



Ml- 
U2* 



-0.2435819 X 10^ 
0.1488425 X lO'^ cm 



cm ^ 
1X6 



Rm 

O-O 

0.1 

0.2 

03 

04 

05 

ae 

07 



Rjnn < J? < Rq 



11.00 A 



-172, 

0.355691862122135882 

0.910756126766199941 
-0.460619207631179620 

0.910227086296958532 
-0.296064051187991117 
-0.496106499110302684 

0.147539144920038962 
-0.819923776793683828 



=u:w 

5.149085 A 
90517 cm~i 
X 10^ 
X 10^ 
X 10^ 
X 10^ 
X 10-* 
X 10* 
X 10^ 
X 10* 



cm 
cm 
cm 
cm 
cm 
cm 
cm 
cm 



-1 
-1 
-1 
-1 
-1 
-1 



^out < ^ 


c/00 


0.0 cm-^ 


Cs 


0.75186131x10^ cm-^A^ 


Cs 


0.1686430x10'' cm^^A** 


Clo 


0.3081961 xlO^o cm-^Ai" 


Aex 


0.40485835 xlO^ cm-^A-T 


7 


4.59105 


/3 


2.36594 A-i 


Derived constants: 


equilibrium distance: 


iJJ= 5.1431(10) A 


electronic term energy 


T;^=-D^= -172.909(40) cm^^ 



ing to the average temperature in the ensemble prepared 
during the experiment. This is not necessarily the cor- 
rect position, because the observed loss of atoms is deter- 
mined by three-body losses. We believe, however, that 
the difference between the resonance position defined by 
the experiment and that of the elastic resonance is not 
more than the experimental accuracy of < 0.01 mT. 

The final potential parameters are given in Table lllll 
andllVlfor the X^T,'^ and a'^Sj states, respectively. The 
large number of digits is only given in order to create 
potentials with the numerical accuracy as applied in the 
fitting process. It does not reflect the physical accuracy 
of the potential. Most of the resonance positions are de- 
scribed well within the experimental uncertainty, see Ta- 
ble HI and thus this approach yields a significant improve- 
ment compared to the ABM model. Additionally, the 
photoassociation data are reproduced with a standard de- 
viation of cr = 0.82 using the measurement uncertainties 
for weights. The fit shows that the i?-dependence of the 
hyperfine coupling is given by Cf — —0.029, Rq = 11.0 
a.u. and AR =1.0 a.u. , where the amplitude c/ is di- 
mcnsionless and we estimate its uncertainty as 30%. Rq 
is chosen such that the variation of the hyperfine cou- 
pling is mainly around the minimum of the triplet state, 
but certainly this choice is a bit arbitrary and different 
choices will alter the amplitude c/. Nevertheless, we find 
that the hyperfine structure decreases from the atomic 
value by a few percent for R < Rq. Similarly, we var- 
ied the parameters for the effective spin-spin interaction 
and obtained aso — —1-56 a.u. and b = 1.0 a.u. , where 



V =14 




100 150 

S(mT) 

FIG. 3: (color online) Magnetic field dependence of the least 
bound I — singlet and triplet states with M = 2. The 
dashed black line represents the dissociation threshold of 
two atoms in the |1, 1) state, whereas the dotted black lines 
represent thresholds for asymptotes F — 1 + F — 2 and 
F = 2 + F = 2 with M = 2. The inset zooms into the mag- 
netic field region around the s-wave resonances at 85.1 mT 
and 90.5 mT. 



b is chosen as in Ref. [il| and the value of 050 has an 
estimated uncertainty of 100%. It is almost a factor two 
smaller than in the above reference, but we should note 
the trend that the second order spin-orbit interaction 
slightly reduces the spin-spin interaction of the atomic 
pair in the molecule. 

Jones et al. psj reported direct determination 
of the dissociation energy Dq of Na2 and found 
i:io/(M=5942. 6880(39) cm-i. We included this value in 
the fitting process and the model potential of X^Y^'^ re- 
produces this energy by 5942.6913 cm~^, thus in agree- 
ment with observation. The dissociation energy of the 
triplet state a'^E+ could not be improved, because it is 
only determined by the less precise spectroscopic data of 
Ref. [4l|. The derived values for Dp, which is the well of 
the Born-Oppenhcimer potential and differs from Dq by 
the zero-point energy, are given in Tablc lTTTl and IIVI 

In Fig.Othe resulting molecular spectrum near the dis- 
sociation threshold is shown for / = states. This spec- 
trum can be compared with that of Fig. [2] to observe the 
infiucnce of the singlet levels via the singlet-triplet cou- 
pling. The most striking feature is the broad avoided 
crossing between the Va = 14 (M5 = —1) and vx = 64 
{Ms = 0) levels around 170 mT, which was already inves- 
tigated by Laue et al. [l^, and used in the ABM fitting 
procedure (see Sec. lIIIB"]) . The effect of the vx = 65 
level is more subtle. The near-degenerate vx = 65 and 
Va = 15 levels are strongly mixed, leading to a signifi- 
cant singlet character of the Va = 15 triplet level. As a 
consequence, the Va = 15 Ms = — 1 state can couple to 
the Va = 14 Ms — 1 state, leading to avoided crossings 
at the intersects of these triplet states. This is high- 
lighted in the inset of Fig. [31 showing the region around 
the 85.1 mT and 90.5mT s-wave resonances. The first 
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<f> -1.86 -0.31 0.42 1.55 0.45 2.94 

' -2 -10 12 3 




50.5 50.6 50.7 50.8 50.9 51.0 51.1 



e(mT) 

FIG. 4: (color online) Coupled-channels results for the group 
oil = A bound states leading to the three observed 5- wave res- 
onances (green triangles), with (green solid lines) and without 
(grey squares) coupling between the different m/ states. For 
the coupled states the expectation value of the operator is 
given on top of the graph, explicitly showing that ruf is not 
a good quantum number for these states. Energy is defined 
relative to the jl, 1) dissociation threshold. 



avoided crossing is narrovif and the effect on the position 
of the resonance caused by the / = 4 state is negligible. 
Hovifever, the second avoided crossing is much broader 
and shifts the position of the resonance caused by the 
/ = 2 state by about 0.6 mT. 

The widths of Feshbach resonances are determined by 
the coupling of the bound levels to the continuum and 
thus directly calculable within the coupled-channels ap- 
proach. TablcH] also lists the widths A defined according 
to Eq.[T]for clastic collisions with incoming s-wave. Only 
two of them, at 90.5 mT and 205.4 mT, are wide enough 
for use in an experiment with magnetic tunability. The 
others are so narrow that the hold time had to be large 
for detecting any enhanced atom loss. In fact, only be- 
cause the applied magnetic field is fluctuating across the 
resonance, do these resonances become observable. In 
this situation the actual resonance width is reflected in 
the experimentally required hold time to observe the loss 
feature, and indeed we notice a clear correlation. We did 
search for the / = 2 d-wave resonance that crosses the 
atomic threshold around 80 mT (see Fig. [2]), but even a 
hold time of 10 s did not result in a loss feature. Accord- 
ing to the coupled-channels calculation the width of this 
particular resonance is even a factor of five smaller than 
the narrowest observed d-wave resonance at 71.6 mT. 

Finally, we focus on the group of ^ = 4 / = 3 [Ms = 0) 
states that causes the observed g-wave resonances. As 
stated earlier, for these states m/ (or mi) is not a good 
quantum number because of strong mixing between the 
different states. The molecular spectrum near the dis- 
sociation threshold is shown in Fig.Ul in which also the 
expectation value {fz) is given. For comparison, the re- 



sults of a calculation that excludes coupling between m f 
states are shown. Both the order of and spacing between 
the states are affected by the coupling. 



IV. CONCLUSION 

We have observed 14 Feshbach resonances for ultracold 
Na in the absolute hyperfine ground state. The improve- 
ment on the field positions of the resonances, now bet- 
ter than 0.04 mT, is very significant. Additionally, the 
newly found g-wave resonances extends the knowledge 
of the asymptotic boimd structure of Na2. A coupled- 
channels calculation describes most of the resonance po- 
sition to well within the experimental uncertainty. The 
i?-dependent hyperfine interaction has a noticeable ef- 
fect on the magnetic field location of the Feshbach reso- 
nance. If we use the atomic hyperfine value for all inter- 
nuclear separations we find that the resonances shifts up 
to 0.05 mT, which is larger than our experimental uncer- 
tainty. With the new model potentials we derive for the 
location of the s-wave resonance in the |1, —1) -I- |1, —1) 
entrance channel 120.75(0.5) mT, in agreement with the 
experimental value of 119.5(2.0) mT 

Our Feshbach spectroscopy data result in a precision 
of about 0.3 MHz in the frequency scale for the weakly 
bound states closest to the atomic asymptote. Thus they 
complement the photoassociation data from the experi- 
ment by Fatemi et al. pj| . which have a precision of 
about 15 MHz. Both data sets are strongly related to the 
triplet state, because the singlet-triplet mixing in Na2 is 
weak as can be seen from the accuracy with which the 
Moerdijk model can describe the position of the Fesh- 
bach resonances. The combined evaluation of all existing 
data gives an accurate description of the a'^Sj potential 
down to about 820 GHz where the measurements of the 
bound states by Ref. [ll| stop. Further down to its min- 
imum only less precise and sparse data exist from the 
spectroscopy by Li Li et al. [4l|. 

TablclVl gives an overview of the most weakly bound 
level energies for the Moerdijk, ABM and coupled- 
channels models. In general there is fair agreement be- 
tween ABM and the coupled-channels approach. One 
finds significant deviations for levels that are only weakly 
linked to the observed Feshbach resonances. 

Ultracold collisions can be modeled with good reliabil- 
ity using the new potentials and the molecular hyperfine 
and second-order spin-orbit coupling. For B = Q the scat- 
tering length for the lowest open channel |l,l) + |l,l)is 
calculated to be 54.54(20)ao (ao = 0.529177 x 10~i°m), 
which is in a good agreement with 55.1(16)ao, derived 
by Crubellier et al. [i^l, but less well to the value 
52.98(40)ao obtained by Samuehs et al. from a spec- 
troscopic study of collision resonances in a molecular 
beam. Earher reported values (see e.g. Refs. [1^, |47| ) 
are less precise but agree with the new value within their 
reported error limits. The small uncertainty for the scat- 
tering length results from the improved precision in the 
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TABLE V: Level energies £^''/h in MHz obtained from 
the different models for the two last bound states of the 
X^T,^ and a'^E^ potentials. The energy reference is the 
atomic pair asymptote |1,1)-|-|1,1). The label * in the column 
ABM indicates that the value was estimated from other spec- 
troscopic sources and not only from Feshbach spectroscopy, 
while the label f denotes that more precise energy differences 
of rotational levels were used (see Sec lIIIBjl . 



X/a V I 


Moerdijk 


ABM 


CC 


a 14 
a 14 2 
a 14 4 


-4976 
-3679 
-764 


-4991(1) 
-3682(2) 
-777(2) 


-4991.5(1.0) 
-3681.4(1.0) 
-779.8(1.0) 


a 15 
a 15 2 




+2014(10)* 
+2479*t 


+2033(15) 
+2358(20) 


X 64 
X 64 2 
X 64 4 




-11000* 
-9444*t 
-5936*t 


-10844(30) 
-9283(30) 
-5749(40) 


X 65 
X 65 2 




-H400* 
-hl904*t 


+1171(30) 
+1825(40) 



surprising given the strong correlation between this scat- 
tering lengths and the one for |1, 1) + |1, 1). The scat- 
tering length of the singlet X^Y.'^ state is mainly de- 
termined by the binding eiiergies of the two last bound 
states vx = 64 and 65 [§, |44| | Thus we could not im- 
prove the precision of the singlet scattering length, but 
its value has changed slightly to 18.81 (80)ao. This fact 
is understandable, because the scattering lengths of the 
singlet and triplet states are correlated by noting that 
the Vx = 65 bound level is significantly mixed with the 
Va ~ 15 one. 

Finally, we propose that the singlet scattering length 
can be improved by looking for a narrow p-wave reso- 
nance in collisions between Na atoms in the |1,1) and 
|1, 0) state. We predict a p-wave resonance at a magnetic 
field of 23.5 mT. It is desirable to improve the molecular 
beam data from about 30 MHz to better than 1 MHz, 
which now seems feasible given femtosecond comb-based 
optical frequency measurements. 



Feshbach resonance locations compared to Refs. [Zj. 
Note that this scattering length is equal to that of the 
|1,-1) + |1, -1) and |1,±1) + |1,0) channels. The other 
scattering lengths relevant for ultracold Na prepared in 
the F = 1 manifold are 52.66(40)ao and 50.78(40)ao for 
the |1, 0) + 11, 0) and |1, 1) + 11, —1) channels, respectively 

We can provide a very reliable value for the scattering 
length of the triplet a'^S^ state, narnely 64.30(40)ao. It 
agrees with the value found in Ref. [43] , but has a 1 ao 
deviation with the value reported in Ref. . This is not 
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